Resolving the Discrepancy of 135 MeV pd Elastic Scattering Cross Sections and 

Relativistic Effects 
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Three precise measurements for elastic pd scattering at 135 MeV/A have been performed with 
the three different experimental setups. The cross sections are described well by the theoretical 
predictions based on modern nucleon-nucleon forces combined with three nucleon forces. Relativistic 
Faddeev calculations show that relativistic effects are restricted to backward angles. This result 
supports the two measurements recently reported by RIKEN and contradicts the KVI data. 

PACS numbers: 21.30.-x, 21.45. +v, 24.10.-i, 24.70.-|-s 



A hot topic of present day few-nucleon system studies 
is to explore the properties of three-nucleon forces (3NFs) 
acting in systems with more than A— 2 nucleons. These 
forces appear the first time in the three-nucleon (3N) 
system where they provide an additional contribution to 
a predominantly pairwise potential energy of three nu- 
cleons. Meson-exchange picture can undoubtedly lead to 
such forces, however they are relatively weak compared 
to nucleon-nucleon (NN) forces, and therefore it is hard 
to approach and also to find evidences for them experi- 
mentally. The first evidence of the 3NF is found in the 
3N bound systems, '^H and ^He. The binding energies of 
these nuclei are not described by exact solutions of the 
three-nucleon Faddeev equations employing modern NN 
forces, e.g. AV18 1|, CDBonn Z\, Nijmegen I, II and 
93 13 (see e.g. Ref. j^). The discrepancy between data 
and theory is explained by adding 3NF, mostly based on 
27r-exchange between three nucleons with the A-isobar 
excitation, such as the Tucson-Melbourne (TM) 5] and 
the Urbana IX 3NF Q. The binding energies show the 
significant contributions of the 3NF, however they could 
only constrain the overall strength. In order to study the 
momentum and/or spin dependence of the 3NF, the 3N 
scattering is one attractive approach; others are to study 
the spectra of nuclear systems up to A=10 0. 

Indication of 3NF for the 3N scattering was first 
pointed out in the cross section minima for nucleon- 
deuteron (Nd) elastic scattering at energies of the in- 
coming nucleon above « 60 MeV by Witala et al. in 
1998 8]. Since then experimental studies of higher- 
energy proton-deuteron (pd) elastic scattering covering 



incident energies of up to w 250 MeV have been per- 
formed extensively and provided precise data of cross 
sections [l3lil|l|JllL4| and spin observables, such as 
analyzing powers [13, IOj ll^ j spin correlation coeffi- 
cients [3 ! and polarization transfer coefficients 0, 0| . 
Precise cross section data for the elastic dp scattering 
taken at RIKEN with 135 MeV/ A deuterons have shown 
large disagreement between data and rigorous Faddeev 
calculations with modern NN forces [ifl, fllL ll^ . Combi- 
nation of these NN forces and 3NFs such as the TM99 [l3| 
(TM99 is a version of the TM force which is more con- 
sistent with chiral symmetry 0, and the Urbana 
IX removed this discrepancy and led to a good descrip- 
tion of the measured cross sections. This result can be 
taken as a clear signature of the 3NF effects in Nd elas- 
tic scattering. However spin observables are not always 
explained by addition of the 3NFs, showing the defects 
in spin parts of the 3NFs 9]. 

The recent measurement of the elastic pd scattering 
with a 135 MeV proton beam and a mixed solid CD2- 
CH2 target at KVI Groningen provided cross sections 
which were in disagreement with the data at RIKEN. The 
KVI data were larger than the RIKEN data by about 10- 
40% and also differed in shape. IF the KVI data were 
correct, the presently available 27r exchange 3NFs would 
be insufficient to explain the difference. Therefore one 
should look for other sources which have not been con- 
sidered to fill the difference, such as 3NFs other than 
27r exchange types and/or relativistic effects or some- 
thing completely new. According to the recent theoreti- 
cal predictions in the framework of the coupled channel 



2 



approach j23|, the ir-p and p-p type 3NFs cause httle 
effect. Therefore relativistic effects could be a candidate 
to fill the difference. We estimate for the first time the 
magnitude of relativistic effects in the 3N scattering. Of 
highest importance, however, is to clarify by experiment 
which are correct, the RIKEN data or the KVI data. 

Aiming to resolve the discrepancy by experiment, we 
performed the following three measurements for elastic 
pd scattering. First, we made a measurement at RIKEN 
with the proton beam and a CD2-CH2 sandwiched solid 
target at the angles where the pp and pd elastic scat- 
tering were simultaneously measured with the magnetic 
spectrograph SMART. Using the well-known elastic pp 
cross sections we can estimate the overall systematic un- 
certainty for the pd cross section. We used the ions of 
270 MeV as the 135 MeV proton beam for convenience of 
acceleration. Secondly, to confirm the angular distribu- 
tion we measured again with 135 MeV/A deuterons. This 
measurement was performed just after the previous pp 
scattering experiment with the same experimental setup 
in order to minimize the systematic uncertainties. Note 
the mass of is almost identical to that of deuteron 
so that we did not need to change any parameters of 
the accelerators or beam transport system. We tried to 
check the fluctuations of the target thickness during the 
experiment by measuring the dp scattering at the fixed 
angle 9c.m. = 69.7°, where the scattered deuterons and 
recoil protons were detected in coincidence in the scat- 
tering chamber. For the same purpose, the cross section 
at dc.m. = 165.1° was measured with the SMART sys- 
tem over several times during the experiment. We also 
measured the carbon background events which were not 
obtained in the previous measurement [ill Il2| . Lastly, 
we performed a totally independent measurement at the 
Research Center for Nuclear Physics (RCNP) of Osaka 
University, using a 135 MeV proton beam and deuterated 
polyethylene target. The absolute normalization of the 
cross sections has been performed by taking data with a 
D2 gas target and the double slit system for which the 
RCNP group has already established the procedure to ob- 
tain the absolute pd cross section |l3j . In this paper we 
would like to present these new data and compare them 
to the old ones and to theoretical predictions including 
relativistic effects to clarify firmly the 3NF effects. 

The measurement of the cross sections with 270 MeV 
beam {p on ^H) was carried out at the RIKEN ac- 
celerator research facility with the SMART system. The 
target was a sandwich of the polyethylene (CII2) with 
thickness of 18.7 mg/cm^ and the self-supporting 99% 
isotopically enriched deuterated poly ethylene foil (CD2) 
with thickness of 21mg/cm^ [Tsl l2l||. The relative devi- 
ation of the CD2 target thickness was estimated to be 
within about 2.5% which was attributed to the inhomo- 
geneity of the CD2 foil. At the lab. angles 10°-14° which 
correspond to the cm. angles 16°-20° for the elastic pd 
scattering the momentum difference of the scattered pro- 
tons from the pp and pd elastic scattering is within 4%. 
In this angular range the protons emitted from each of 



these reactions were measured simultaneously with the 
SMART spectrograph and their energy spectra were com- 
pletely resolved. The yields for the elastic pd scattering 
were obtained by subtracting carbon contributions in the 
excitation energy spectra. The elastic pp scattering yields 
were obtained by subtracting the backgrounds from the 
p-\-d^p-\-p-\-n breakup reaction and from the proton 
scattering on the carbon. For interpolation purposes the 
breakup background contribution was assumed to be a 
third-order polynomial in the proton energy. The mea- 
sured cross sections for pp elastic scattering were com- 
pared with the values calculated by the phase-shift anal- 
ysis code SAID Q2i and found to be consistent within 2%. 
Thus we estimated the overall systematic uncertainty of 
the measured cross section data to be 3% at most. 

We made the cross section measurement with 270 
MeV deuteron beam (d on ^H) in the angular range 
6»e.m. = 10°-180°. The CH2 target used in the proton 
beam experiment was employed as a hydrogen target. 
For the forward scattering [Ocm. < 90°) the scattered 
deuterons were detected while for the backward scatter- 
ing {Oc.ra. > 90°) the recoil protons were measured. The 
statistical errors of the cross sections are within 1.6%. 
The fluctuation of the target thickness is within 3%. The 
uncertainty due to carbon background subtraction is less 
than 5%. The overall systematic uncertainties which in- 
clude also the uncertainties by pp experiment are esti- 
mated to be 6%. 

The experiment performed at RCNP used a proton 
beam in conjunction with the high resolution spectrome- 
ter Grand Raiden. The proton beam was accelerated up 
to 135 MeV by the AVE and ring cyclotrons and bom- 
barded the same CD2 foil we used in the experiment at 
RIKEN. The proton beam was stopped in a Faraday cup 
in the scattering chamber, except for the gaseous target 
measurement at the angle ^lab. = 25.5°. In this case 
the beam was stopped in a Faraday cup located down- 
stream outside the scattering chamber. The scattered 
protons or deuterons were momentum analyzed by the 
Grand Raiden. The protons were measured at the an- 
gles 6'c.in. < 90° and the recoil deuterons were detected 
at the angles 6c.m. > 90°. The measured angles were 
6ie.„i. = 17.0°-157.7°. The yields from D2 were obtained 
by subtracting carbon contributions in the excitation en- 
ergy spectra. To normalize cross sections taken with the 
CD2 target a measurement with a D2 gas target was per- 
formed at the lab. angles 25.5° and 60° which corre- 
sponded to the cm. angles 39.0° and 87.4°, respectively. 
The D2 gas target was contained in the cell of a cylinder 
of 40 mm diameter made of 200 /zm-thick aluminum. The 
absolute gas pressure was continuously monitored with a 
precision better than 0.2%. The temperature of the tar- 
get cell was checked during the measurement and kept 
at room temperature. A double-slit system was used to 
determine precisely the target volume and the solid angle 
of the Grand Raiden spectrometer. The effective target 
thickness and the solid angle were calculated by Monte 
Carlo simulations. Spectra with the empty gaseous cell 
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were also measured to determine background contribu- 
tions from the aluminum cell. An additional measure- 
ment was performed with hydrogen gas replacing the 
deuteron one in order to cross-check the experimental 
setup at the angle 6iah. — 25.5°. The measured cross 
section of pp scattering is consistent within 3% with the 
calculated results by the SAID program. The statisti- 
cal errors of the pd elastic scattering cross sections are 
smaller than 1.4%. The absolute normalization was es- 
timated to be 3% by elastic pp scattering measurement. 
The uncertainty due to carbon background subtraction 
for the excitation energy spectrum is 3%. There is also 
the uncertainty of 2.5% attributed to the inhomogeneity 
of the CD2 foil. The overall systematic uncertainties are 
estimated to be 5% at most. 

All the experimental results are shown in Fig. ^ The 
data taken at RCNP are shown with solid diamonds. 
The open squares (circles) are the data measured with 
the proton (deuteron ) b eam at RIKEN. The data pub- 
lished in Refs. [iflli3jli3 shown with open diamonds. 
The KVI data reported in Ref. |Q| are shown with open 
triangles. Only statistical errors are presented. The 
very good agreement between the independent measure- 
ments allows us to conclude that the systematic uncer- 
tainty due to the detection setup is small. Comparison 
of the new three sets of data to the previously reported 
ones 0,0,01 supports the previous measurements and 
shows a clear disagreement with the KVI data [T^ . 

In Fig.n]we compare the data to theoretical predictions 
based on modern NN forces and their combinations with 
3NFs m m m we used the modern NN forces 
AV18, CDBonn, Nijmegen I and II combined with the 
3NF TM99 with the cut-off A values which lead for 
a particular NN force combined with the 3NF TM99 to 
a reproduction of the binding energy. In case of the 
AV18 NN force we also combined it with the 3NF Urbana 
IX . 

As is seen in Fig. ^ various NN force predictions (light 
shaded band) clearly underestimate the cross section 
data for the angles dc.m. > 90°. The narrow band of 
predictions reflects the weak dependence on the choice of 
the nearly on-shell equivalent NN interaction. The inclu- 
sion of the TM99 (the dark shaded band) or, in case of 
AV18, of the 3NF Urbana IX (solid line), leads to a very 
good description of the data. 

In view of such a good description of the cross sec- 
tion it is interesting to find out how large relativistic 
effects are at 135 MeV/A. This was done assuming that 
only NN forces are acting. We followed a formalism for 
treating the relativistic three-body Faddeev equations of 
Ref. l27>| with a boosted two-nucleon potential V ex- 
pressed in terms of the relativistic potential v given in 
the NN cm. system, 



menta of the nucleons in their NN cm. {uj{k) 



V{P) 



p . (1) 



2-\/ fc 2 _|_ 77j2~)^ -y^g j-^Q^ calculate the matrix elements 
of the boosted potential in all its complexity but re- 
stricted only to the leading order terms in a P/w and v /lo 
expansion. 



V{k,k';P) = v{k,k')\ 1 



P^ 



2uj{k)uj{k') 



(2) 



The momentum P is_|;he total ^momentum of the two- 
nucleon system, and k and — k are the individual mo- 



We checked the quality of the approximation of Eq. (|2Jl by 

calculating the deuteron wave function (pdik) when the 
deuteron is moving with a momentum corresponding to 
135 MeV/A. The resulting deuteron binding energy and 
the deuteron Z3-state probability for the deuteron in such 
a motion are close to the values for the deuteron at rest. 
Choosing the exact expression Eq.(l) those properties 
come out exactly. A relativistic potential v was generated 
from the nonrelativistic NN potential AV18 by perform- 
ing the scale transformation of Ref. This should 
be improved in future studies by allowing for dynami- 
cally dictated relativistic features in the potential. Using 
V{ k , k' ; P) the boosted t-matrix elements t{k , k' ; P) 
are calculated and they form the dynamical input for the 
3N Faddeev equation with the relativistic form of the 
free propagator Go [s^]. To solve this equation in the 
relativistic case it is most convenient to use instead of 
the standard Jacobi momenta the relative momen- 
tum k in the NN cm. subsystem and the momentum 
~q of the spectator nuclcon in the 3N cm. system. In 
the nonrelativistic limit the momentum k reduces to the 
standard Jacobi momentum "p^ 2f| . The relativistic for- 
mulation applied is of the Bakamjian-Thomas type and 
belongs to the instant form of relativistic dynamics lU . 

Nowadays partial wave decomposition is still required 
to solve numerically 3N Faddeev equations. The stan- 
dard partial wave states |2^, however, are generalized 

due to the choice of the NN-subsystem momentum k 
and the total spin s both defined in the NN cm. system. 
This lead to Wigner s pin rotations when boosting to the 
3N cm. system |3l|, |32, resulting in a more complex 
form for the permutation matrix element "s^l than used 
in j26| . The details of our relativistic formulation and its 
numerical performance are given in Ref. • 

A restricted relativistic calculation with j < 2 par- 
tial waves states showed that Wigner spin rotations have 
only negligible effects on the cross section at 135 MeV/A. 
Thus when performing the fully converged calculation 
(j < 5, J < 25/2) we neglected the Wigner rotations 
completely. The resulting cross sections are shown in 
Figffl It is seen that the effects of relativity are visi- 
ble only in the backward angular region for 0c. m. > 160° 
where they increase the cross section by up to about 15%. 
For Oc.m. < 160° the effects of relativity are practically 
negligible. 

Summarizing, we performed measurements of the cross 
sections for elastic pd scattering using 135 MeV proton 
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and 270 McV dcutcron beams. Present experimental ar- 
rangements allowed us to get precise cross section data. 
The three sets of data presented here taken with three 
different experimental setups completely support the pre- 
vious measurement and disagree with the KVI data. The 
agreement of our new sets and our old set, measured with 
different experimental setups, gives confidence that the 
systematic errors are small. 

Comparison of our data with theory based on different 
NN forces combined with current 3NFs revealed clear ev- 
idence for the action of 3NFs. The discrepancies between 
the 135 MeV cross section data and the pure two-nucleon 
force predictions can be removed by including the 3NFs 
TM99 or Urbana IX. The conclusion that 3NF effects 
are seen in the region of the cross section minimum is 
further supported by including relativity in the instant 
form of relativistic dynamics as proposed by Bakamjian 
and Thomas. This leads to small relativistic effects at 



backward angles, but negligible contributions in the min- 
imum, leaving 3NFs as the only plausible mechanism to 
resolve the discrepancies between NN theory and data. 
Our results clearly indicate the usefulness of Nd elastic 
scattering for the study of 3NFs. Nd elastic scattering 
cross sections together with spin observables at higher 
energies, where there are still discrepancies, will provide 
an important information to test forthcoming additional 
3N force mechanisms. 
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FIG. 1: (Color online) Differential cross section for elastic Nd 
scattering at 135 MeV/nucleon. The light shaded band con- 
tains NN force predictions (AV18, CDBonn, Nijmcgcn I and 
II). The dark shaded band results when tiiey are combined 
with the TM99 3NF. The solid line is the AV18 + Urbana 
IX prediction. The dashed and dotted lines are the results 
of relativistic and nonrelativistic Faddeev calculations based 
on the NN forces AV18. The symbols are data from different 
measurements (see text). 
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